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Abstract. Using the Darboux method and its relation with supersymmetric quantum 
mechanics we construct all SUSY partners of the harmonic oscillator. With the help of 
the SUSY transformation we introduce ladder operators for these partner Hamiltonians 
^ \ and shown that they close a quadratic algebra. The associated coherent states are 

■ constructed and discussed in some detail. 
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Since the early days of quantum mechanics there has been enormous interest in 
exactly solvable quantum systems. In fact, Schrodinger himself initiated a program 
[1] which resulted in the famous Schrodinger-Infeld-Hull factorization method [2]. 
^ \ In the last 10-15 years this program has been revived in connection with super- 

symmetric (SUSY) quantum mechanics [3]. To be a little more precise, it has been 
O^. found [4] that the so-called property of shape-invariance of a given Schrodinger 

potentials, which is in fact equivalent to the factorization condition, is sufficient 
for the exact solvability of the eigenvalue problem of the associated Schrodinger 



;_j ■ Hamiltonian. However, SUSY quantum mechanics has also been shown to be an 



effective tool in finding new exactly solvable systems. Here in essence one utilizes 
the fact that SUSY quantum mechanics consists of a pair of essentially isospectral 
Hamiltonians whose eigenstates are related by SUSY transformations. This is the 
basic idea of a recent construction method for so-called conditionally exactly solv- 
able potentials [5] . Here one constructs a SUSY quantum system for which, under 
certain conditions imposed on its parameters, one of the SUSY partner Hamiltoni- 
ans reduces to that of an exactly solvable (shape-invariant) one. Other approaches, 
which are also based on the presence of pairs of essentially isospectral Hamiltonians, 
go back to an idea formulated by Darboux [6], are based on the inverse scattering 
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method [7], or on the factorization method [8]. Clearly, these approaches are closely 
connected to each other and to the SUSY approach. 

In this paper wc will construct with the help of the Darboux method all pos- 
sible SUSY partners of the harmonic oscillator Hamiltonian on the real line and 
discuss their algebraic properties in some detail. In doing so we review in the next 
section the Darboux method and explicitly show its equivalence to the supersym- 
metric approach. Section 3 then briefly presents the basic idea for the construction 
of conditionally exactly solvable (CES) potentials. Section 4 is devoted to a de- 
tailed discussion of the harmonic oscillator case. Here we first present all possible 
SUSY partners of the harmonic oscillator and give explicit expressions for the cor- 
responding eigenstates. Secondly, with the help of the standard ladder operators of 
the harmonic oscillator we introduce similar ladder operators for the SUSY part- 
ners and show that they close a quadratic algebra, which is also briefly discussed. 
Finally, we introduce so-called non-linear coherent states which are associated with 
this non-linear algebra. The properties of these coherent states are discussed in 
some detail. 



THE DARBOUX METHOD 

In this section we briefly review the Darboux method [6] and show its connection 
to supersymmetric quantum mechanics [3]. In doing so we start with considering 
a pair of standard Schrodinger Hamiltonians acting on L"^' 



and a linear operator 



h d 

A = ^= — + ^ix) , $ : K ^ R , (2) 

obeying the intertwining relation 

H+A = AH^ . (3) 

It is obvious that this intertwining relation cannot be obeyed for arbitrary functions 
V± and In fact, the relation (3) explicitly reads 



'2m 



^"{x) + V+{x)^{x) - VL{,x) - ^{x)V-{x) \ 1 = 

(- ^'{x) + ^V.{x)-^ V+ix)] I 
V2m V2m J da 



(4) 



As the unit operator 1 and the momentum operator (i.e. d/dx) are linearly in- 
dependent, their coefficients have to vanish. In other words, we are left with two 
conditions between the three functions V± and $: 



V.{x)^V+{x)-^<^>'{x), (5) 
v2m 

<^"{x) + V+{x)<^>{x) - VL{x) - <^>{x)V_{x) = . (6) 

2m V2m 

Inserting the first one into the second one and integrating once we find 

h 



/2m 



^'{x)-V+{x)+^^{x)^-e , (7) 



where e is an arbitrary real integration constant sometimes called factorization 
energy [3] . With this relation and with (5) we can express the two potentials under 
consideration in terms of the function 

V±{x) = ^\x) ± ^'{x) +s. (8) 
v2m 

At this point one realizes that these are so-called SUSY partner potentials [3]. In 
fact, using relations (8) we note that 

H+ = AA^ + e, H_ = A^A + e. (9) 

These supersymmctric partner Hamiltonians are due to the intertwining relation 
(3) essentially isospectral, that is, 

spec H+\{e} = spec H^\{e} . (10) 

Their eigenstates are related via SUSY transformations. To make this more explicit, 
let us denote by the eigenstates of H± for eigenvalues > e, 

^±l0^)=i^n|0^). (11) 

Then these states are related by SUSY transformations [3] 

\<i^:)^^^=A\<p-), i0.) = ^=^^V;t)- (12) 

In addition to the states in (11) one of the two Hamiltonians H± may have an addi- 
tional eigenstate with eigenvalue e obeying the first-order differential equation 
A\(f)~) = and ^4^10+) = 0, respectively. In terms of the function $ they explicitly 
read 

(/)^(x) = A^±exp|±^^ J dx$(x)| , (13) 

where N± stands for a normalization constant. Clearly, only one of the two solutions 

(13) may be square integrablc. This situation corresponds to an unbroken SUSY. 
If none of them is square integrable then SUSY is said to be broken [3] . 



The Darboux method reviewed in this section can now be used to find for a 
given potential, say V+, all its possible SUSY partners Firstly, one has to 

solve equation (7), that is, finding all possible SUSY potentials $. This in fact 
corresponds to find all possible factorizations for the corresponding Hamiltonian 
Finally, the corresponding SUSY partner V- can be obtained via (5). In this 
way one can construct new exactly solvable potentials. The parameters involved 
in the SUSY potential turn out to obey certain conditions and therefore these new 
potentials are more precisely called conditionally exactly solvable (CES) potentials. 
Let us note that the Darboux method may be generalized to intertwining operators 
containing higher orders of the momentum operator [9] . 



In this section we give some more details on the construction of CES potentials 
using the Darboux method. As just mentioned above we start with a given potential 
V+ and try to find all its associated SUSY potentials. That is, we have to find 
the most general solution of the generalized Riccati equation (7). In doing so we 
will first linearize this non-linear differential equation via the substitution = 
{h/ y/2m)u'{x) /u{x), 



which is actually a Schrodinger-hke equation for V^. Note, however, that we are 
not restricted to normahzable solution of (14). In other words, the energy-like 

parameter e is up to now still arbitrary. 
In terms of u the linear operator A reads 



and thus is only a well-defined operator on L^(M) if u does not have any zeros on 
the real line. As a consequence we may admit only those solutions of (14) which 
have no zeros. Form Sturmian theory we know that this is only possible if e is below 
the ground-state energy of which we will denote by Eq. Hence, we obtain a 
first condition on the parameter e, which reads e < Eq. This also implies that e 
does not belong to the spectrum of In fact, the associated eigenfunction (13) 
would read (t>'^{x) — N^u{x), which is not normahzable due to condition put on e. 

The above condition on e is still not sufficient to guarantee a nodeless solution. 
Being a second-order linear differential equation (14) has two linearly independent 
fundamental solutions denoted by Ui and U2- Hence, the most general solution for 
£ < £'o is given by a linear combination of the fundamental ones: 



MODELLING OF CES POTENTIALS 



2m 



u"{x) -\- Vj^{x)u{x) = su{x), 



(14) 




(15) 



u{x) = aui{x) + I3u2{x) 



(16) 



Therefore, the condition that u does not vanish also imposes conditions on the 
parameters a and /3, which have to be studied case by case [5]. 

Let us now assume that is an exactly solvable Hamiltonian, which means 
that its eigenvalues En and eigenstates are exactly known in closed form. For 
simplicity we have assumed that if+ has a purely discrete spectrum enumerated 

by n — 0,1,2,... such that e < Eq < Ei < Then via the method outhned 

above one can construct all its SUSY partners H_ which are conditionally exactly 
solvable due to the conditions which have to be imposed on the parameters a, (3 
and e. By construction the eigenvalues of if+ are also eigenvalues of H_ and the 
corresponding eigenfunctions are obtained via the SUSY transformation (12). In 
the case of unbroken SUSY H_ has one additional eigenvalue e which belongs to 
its ground state given by 4>j{x) — N^/u{x). Finally, we note that in terms of u 
the partner potentials read 

V_{x) = ^('^\ -V^{x) + 2e (17) 
m \ u{x) ) 

and form a two-parameter family label by e and /3/q;. Note that only the quotient 
13 /a ox its inverse is relevant for (17). For various examples of CES potentials 
found by this method see [5]. Here we limit our discussion to those related to the 
harmonic oscillator. 



THE HARMONIC OSCILLATOR 



In this section we will now construct all possible SUSY partner potentials for 
the harmonic oscillator V"+(x) = {m/2)uj^x'^ , u; > 0, via the Darboux method. The 
corresponding Schrodinger-like equation (14) reads in this case^ 

1 1 

-- u"{x) + - x^u{x) = eu{x) (18) 

and has as general solution a linear combination of confluent hypergeometric func- 
tions 

u(x) = e-^y'[a^F,C-^,lx')+Px^F^C--^,l,x')] . (19) 

The condition that u docs not have a real zero implies that a must not vanish and 
thus can be set equal to unity without loss of generality. Furthermore, /3 has to 
obey the inequality [5,10] 

l/3|</3c(£):=2^y^. (20) 



Prom now on we will use dimensionless quantities, that is, x is given in units of ^Jhjmu) and 
all energy- like quantities are given in units of /kj. 



The corresponding partner potentials of the harmonic oscillator then read according 
to (17) 



u\x) 



u[x) 



1 2 ^ 

-x^ + 2e . 
2 



(21) 



We note that for the above u SUSY remains unbroken and therefore, the spectral 
properties of are given by 

spec H_ — {e, Eq, Ei, . . .} with E^ = n + \ , n = 0, 1, 2, . . . , 



3^ 3 



:22i 



exp{-xV2} 



[Vvr2'^+in!(n + l/2-e)]V2 



n+l 



x) + 



u\x) 
u{x) 



- X HJx) 



where if„ denotes the Hermite polynomial of degree n. Figures of the potential 
family (21) for various values of e and {3 can be found in [5]. Here let us stress 
that one can even allow for complex valued (3 e C\[— /3c(£), /3c(^)] which in turn 
will give rise to complex potentials generating the same real spectrum [10]. We also 
note that the present CES potential (21) contains as special cases those previously 
obtain by Abraham and Moses [7] and by Mielnik [8]. See also [5] for a detailed 
discussion. 



Algebraic Structure 

We will now analyse the algebraic structure for the partner Hamiltonians of the 
harmonic oscillator. In fact, using the standard raising and lowering operators of 
the harmonic oscillator H+ — AA^ -\- e — a)a-\-l/2, 

which close the linear algebra 

a] = -a , [//+, o)\ = , [a, a^] = 1 , (24) 

one may introduce via the SUSY transformation (12) similar ladder operators for 
the SUSY partners [11] 

B = A'^aA , = A^a^A , (25) 



which act on the eigenstates of H_ in the following way 



^I0n+i) = \l{n+\-e){n+l){n+l-e)\(t>-) , 



St|0-) = ^{n+\-e){n+l){n+\-e)\ct>-+^) , 
S|0o) = O, S|4-) = 0, St|4-) = 0. 



(26) 



The last two relations explicate that the ground state \4>~) of H_ is isolated in 
the sense that it cannot be reached via B from any of the excited states and, vice 
versa, the excited states cannot be constructed with B'^ from l^j). These ladder 
operators close together with the Hamiltonian i7_ the quadratic, hence non-linear, 
algebra 



S] = -S , B^\ = B^ , [S, B^\ = 3i/! - 4.eH_ + . (27) 



This quadratic algebra belongs to the class of so-called W2 algebras and may be 
viewed as a polynomial deformation of the sm(1, 1) Lie algebra. Such deformations 
have been discussed by Rocck [12] and, within a more general context, by Karassiov 
[13] and Katricl and Quesne [14]. The quadratic Casimir operator associated with 
the algebra (27) reads 

C ^ BB^ -<lf{H_) , ^>{H_) -^>{H_-1) = ZHl-AeH_ + e^ . (28) 

In the Fock space representation (26) we have the following explicit expression 



and the relations BB^ = ^(i/_) and B'^B = ^(i/_ - 1). Hence the Casimir (28) 
vanishes within this representation as expected [13,14]. 



Let us now construct the non-linear coherent states [15] associated with the 
quadratic algebra (27). There are several ways to define such states [16]. Here we 
will define them as eigenstates of the "non-linear" annihilation operator B, leading 
essentially to so-called Barut-Girardello coherent states [17]. We also note that the 
construction procedure presented below is very similar to that of coherent states 
associated with quantum groups [18]. 

Let us note that the ground state of if_ is isolated and therefore we may 
construct the coherent states over the excited states {\(pn)}neNo only. For this 
reason we make the ansatz 



*(i/_) = {H_ - e){H_ + \){H_ + l- e) 



(29) 



Non-linear coherent states 



00 




(30) 



n=0 



where ji is an arbitrary complex number and the real coefficients c„ are to be 
determined from the defining relation 



B\^x)^^x\^x)^Y.'',,^x^B\(t>-) . (3i) 



n=0 



Using relations (26) we obtain the following recurrence relation for the c„'s, 

Cn+i = c„ [(n + 1 - e) (n + 1) (n + I - e)] . (32) 
That is, the coefficients c„ for n > 1 can be expressed in terms of cq, 

Cn = Co [n!(i -£)„(§ -£)„]"'/' (33) 

where {z)n = T{z + n)/r{z) denotes Pochhammcr's symbol. The remaining coeffi- 
cient Co = co(/x) is determined via the normalization of the coherent states 

°^ |,,|2n 1 
n=0 ^2 '^)n\2 

Thus, we can express co in terms of a generalized hyper geometric function [19] 

Co'(/x)= oF2{\-e,l-s-M?) ■ (35) 

Let us now discuss some properties of these non-linear coherent states. First we 
note that these states are not orthogonal lor ji^ v as expected: 

(/i|z/) = Co(/i) Co(z/) (I - e, I - £; //V) . (36) 

Secondly, let us investigate whether these states form an over complete set. In other 
words, we consider the question: Can these states generate a resolution of the unit 
operator? For this we have to recall that the non-hnear coherent states have been 
constructed over the excited states of Therefore, we start with postulating 
a positive measure p on the complex //-plane obeying the following resolution of 
unity: 

/ dp{ii\ii)\ii){ii\^i-\ct>:){(t>:\. (37) 

Within the polar decomposition p — y/x e"'^ we make the ansatz 

with a yet unknown positive density a on the positive half-line. Inserting this 
ansatz into (37) we obtain the following conditions on a 



f 

Jo 



.V r(i -e + n)r(| -e + n) 
dxa{x)x^^r{n + l) ^^p(i_^jp|l_^^ ' , n = 0,l,2,.... (39) 
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FIGURE 1. The radial density f{x) = a{x) / (^{^/x) giving rise to the resolution of unity (37) 
with (38) as a function of a; = and for various parameters £ < \- 

Hence, o" is a probability density on the positive half-line defined by its moments 
given on the right-hand side of (39). Let us note that the integral in (39) may be 
viewed as a Mellin transformation [20] of a and in turn the latter is given by the 
inverse Mellin transformation of the moments. This inverse Mellin transformation 
turns out to lead to the integral representation of Meijer's G- function [19]. In other 
words, we have the explicit form: 

= r(i_,)r(|-£) (^|0' -i - I - ^) • (40) 

In Figure 1 a plot of the radial density = 2t: dp{^* , ^) / {dipd\^\'^) is given 

showing that it leads to a well-behaved positive measure on the complex ^-plane. 

Finally, let us point out that similar non-linear coherent states associated with the 
CES potentials of the radial harmonic oscillator have been constructed in [15]. In 
that case broken as well as unbroken SUSY can be considered and the corresponding 
symmetry algebra is a cubic one. In analogy to the discussion in [15] one can show 
that the coherent states discussed here are also minimum uncertainty states. 
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